/M-/// 0/ 




oh 1 




A simulation analysis of the fate of phytoplankton 
within the Mid-Atlantic Bight 



{NASA-CB- 177265) A SIMULATION ANALYSIS OF N86-26670 

THE FATE OF P BY TOPI A N K10N WITHIN THE 
MID-ATLANTIC EIGHT {University of South 

Florida, St. Petersburg.) 62 p Unclas 

HC A04/MF A0 1 CSCL 08A G3/43 43604 


by 


John J. Walsh, Dwight A. Dieter le, and Mark B. Meyers 


Department of Marine Science 
University of South Florida • 
St. Petersburg, Florida 33701 




ABSTRACT 


A time-dependent, three-dimensional simulation model of 

wind-induced changes of the circulation field, of light and nutrient 

regulation of photosynthesis, of vertical mixing as well as algal 

sinking, and of herbivore grazing stress, is used to analyze the 

seasonal production, consumption, and transport of the spring bloom 

within the mid-Atlantic Bight. Our case (c) of a 58 -day period in 

February-April 1979, simulated primary production, based on both 
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nitrate and recycled nitrogen, with a mean of 0.62 g C m day over 

the whole model domain, and an export at the shelf-break off Long 
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Island of 2.60 g chi m day within the lower third of the water 
column. About 57% of the carbon fixation was removed by herbivores, 
with 21% lost as export, either downshelf or offshore to slope waters, 
after the first 58 days of the spring bloom. Extension of the model 
for another 22 days of case (c) increased the mean export to 27%, 
while variation of the model's parameters in 8 other cases led to a . 
range in export from 8% to 38% of the average primary production. 
Spatial and temporal variations of the simulated algal biomass, left 
behind in the shelf water column, reproduced chlorophyll fields sensed 
by satellite, shipboard, and in situ instruments. 



INTRODUCTION 


Satellite time series of the 1979 spring bloom within the 

mid-Atlantic Bight (WALSH, DIETERLE and ESAIAS, 1987a) and moored 

fluorometer time series of the 1984 spring bloom in the same region 

(WALSH, WIRICK, PIETRAFESA, WHITLEDGE, HOGE and SWIFT, 1987b) 

suggested that the daily export of algal biomass to slope waters might 
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be 1.8-2. 7 g chi m day across the shelf -break. Such time series 
lacked depth resolution, however, since the CZCS sampled, at most, the 
upper 10 m of the water column and the in situ instruments sampled 
perhaps the lower 5 m at the 80 m and 120 m isobaths. Furthermore, 
the algal biomass detected by the CZCS over space and by the moored 
fluorometer over time only represented what was left behind in the 
water column by the bacteria, zooplankton, and benthos of this shelf 
ecosystem. 

To place these daily estimates of export in the context of 
seasonal production, consumption, and transport of the spring bloom, a 
simulation model of the plankton dynamics, over 6000 grid points in 
the horizontal plane and 3 depth layers in the vertical, was run for 
58 days from 28 February to 27 April 1979 in response to wind-induced 
changes of the barotropic circulation field. The details of the 
time-dependent terms for the physical circulation, for the light and 
nutrient regulation of photosynthesis, for the vertical mixing, for 
the algal sinking, and for the grazing stress at -18000 grid points 
are discussed below. 

During the 6 yr interval between these CZCS and fluorometer time 
series, a set of biochemical surveys was also conducted each April 
aboard the R/V Kelez , Advance II , and Cape Henlopen . The April 



CHLOROPHYLL 


Figure 1. The mean chlorophyll content (mg m' 3 ) of surface waters 
(0-30 m) of the mid-Atlantic Bight at 27 stations 
re-occupied each April from 1980 to 1984 (W. GREGG, 

personal communication). 


1980-84 chlorophyll content over the first attenuation depth (37% 
light level) at the same 27 station locations in estuarine ( <20 m), 
shelf (20-200 m), and slope (>200 m) regions of the mid-Atlantic Bight 
is plotted (W. GREGG, personal communication) as Figure 1. The mean 
chlorophyll content of these waters may have doubled over 5 yr, from 
2.38 pg chi l * in April 1980 to 5.25 pg chi l * in April 1984. Not 
surprisingly, the chlorophyll content of slope waters was the same 
each April, with the possible increase of algal biomass taking place 
in shelf and estuarine waters. 

In contrast to pristine river contents of -5-10 pg-at N ^ 
(WALSH, 1984), the dissolved nitrogen content at the freshwater end of 
the Hudson, Delaware, and Susquehanna estuaries was >50 pg-at N l ^ in 
1980 and may still be increasing as a result of sewage wastes and 
agricultural leachate. The volume of these estuaries, and thus the 
residence time for algal uptake of dissolved nitrogen, increases 
southward, with consequent export to the adjacent shelf of mainly 
dissolved nitrogen from the Hudson River plume and of mainly 
particulate nitrogen from the Chesapeake Bay plume (Fig. 2). 

Most of the chlorophyll in Chesapeake Bay occurs upstream of the 
25X isopleth of salinity, while most of the algal biomass in the 
Hudson River system is produced downstream of this isopleth (Fig. 2). 
At the 20 m isobath, seaward of the Hudson River estuary, tenfold more 
nitrate is thus present during the spring bloom than at the same depth 
off Delaware or Chesapeake Bay (Fig. 3). An assessment of the 
estuarine contribution to algal export at the shelf-break was also 
examined in the simulation model by altering the boundary flux of 
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Figure 2. The nitrate and chlorophyll distributions within the 
a) Hudson, b) Delaware, and c) Chesapeake estuaries as a 
percentage function of their maximum value and salinity 
during winter (o) and summer (x) - after CARPENTER et al. , 
1969; McCarthy et al., 1977; DECK, 1981; and SHARP et al., 
1982. 








nutrients at the mouths of these three estuaries as well as that from 
Long Island Sound (Fig. 4). 


METHODS 

1 . Circulation 

Ignoring tidal forces and the horizontal stress terms of the 
Navier-Stokes equations, their cross-differentiation and the 
continuity equation leads to a steady-state, depth- integrated form of 
the vorticity equation as 
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where the first term is a description of vortex shrinking or 
stretching as a water parcel crosses an isobath, in which H is the 
bottom depth, and <|> is the sea level potential, defined as the product 
of the acceleration of gravity and the height of the free surface 
elevation. The second and third terms are respectively the curl of 
the components of the bottom stress, and B x » and of the wind 
forcing, F^ and F x in a left-hand cartesian coordinate system, with z 
positive downwards. 

The bottom stress can be defined (HOPKINS and DIETERLE, 1983) in 
terms of the sea level potential and the wind stress by 
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Figure 3. The seasonal structure of nitrate (yg-at SI' 1 ) at the 20 m 
isobath off the a) Hudson, b) Delaware, and c) Chesapeake 
estuaries during March-September (T. WHITLEDGE, personal 
communication). 
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where and b 2 represent that portion of the bottom stress due to the 
wind-driven velocity components. In waters deeper than the surface 
Ekman layer (H £ 30 m), the bottom stress terms become similar to 
those of HSUEH and PENG (1978), e.g.. 


B y « a ( cos 8) .|£+ a (sin 8) (4) 

where 0 is a bottom veer angle, and a represents a scaling depth. It 
is somewhat analogous to the depth of the bottom Ekman layer, in which 
a might be 5-15 m, with 0 ranging from 0° to 45° (HSUEH and PENG, 

1978; HSUEH, 1980; HAN, HANSEN and GALT, 1980; and CSANADY , 1981). 

Substituting eqs. (2) and (3) in eq. (1), assuming that the wind 
curl is negligible, and setting h = H + a^, leads to 
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which can be solved as a steady-state, boundary-value problem for the 
sea level potential. 

Wind forcing also enters the solution of eq. (5), at the coastal 
boundary condition, where the depth- integrated transport normal to the 
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Figure 4. A 6,000 element curvilinear grid of a coupled physical- 
biochemical model of the mid-Atlantic Bight in relation to 
upstream (N1-N6) and interior (LI 1 , L13, LTM, N23, N22, 
N31, N41, and N32) current meter moorings. 



coastline, U = / ■ udz, is assumed to be zero, except at the mouths of 

H 

the three estuaries. The alongshore transport, V = / v dz, at the 
coastal boundary (Fig. 4) is governed by the depth- integrated shallow 
water equation 
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where all the terms are defined above. At steady-state, i.e., — = 0, 

O t 

eq. (6) becomes 
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except at the estuaries, where U is non-zero. Use of eq. (3) results 
in 
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at x = 0. For 13 time intervals between 28 February and 27 April 
1979, the alongshore component of the mean wind stress observed at 
John F. Kennedy airport (Fig. 5) was entered in the simulation, with a 
new steady-state solution of ()) calculated (HOPKINS and DIETERLE, 1983) 
at each grid point (Fig. 4) from a finite difference form of eqs. (5) 
and (8). 

The open water boundary conditions for must also be specified, 
of course, for solution of this boundary value problem. The sea 
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Figure 5. The daily and mean wind stress (dynes cm" 2 ) over 13 
periods as measured at the John F. Kennedy airport between 
25 February and 1 May 1979. 


elevation at the upstream boundary of the grid was computed from the 
observed currents (BEARDSLEY, CHAPMAN, BRINK, RAMP and SCHLITZ, 1985) 
at N1 to N5 (Fig. 4) for each time interval (Fig. 5). Since the April 
1979 currents were mainly barotropic in the mid-Atlantic Bight 
(HOPKINS and DIETERLE, 1987) and hydrographic data on density 
gradients were sparse, we did not attempt to subtract the baroclinic 
component of flow from the current meter records. At the offshore 
boundary, sea elevation was set to a constant 5 cm, i.e., mass 
exchange was confined at the shelf-break to the surface and bottom 

9 2 d) 

Ekman layers. At the downstream boundary (Fig. 4), — f = 0 was 

dy 

specified (HSUEH, 1980). 

At each point of the curvilinear grid (HOPKINS and DIETERLE, 

1983) of the model, i.e., a length of 886 km alongshore at the coast 
and 676 km at the shelf -break (Fig. 4), <J> was first computed for each 
of the 13 wind cases and then inserted in the steady- state Ekman 
equations 
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for evaluation of u and v. Analytical solution of eqs. (9) and (10) 
then allowed (HOPKINS and SLATEST, .1986) continuous depth resolution 
of the flow field for different values of the wind stress and vertical 
eddy diffusivity, K z , at each steady state. In the biochemical 
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calculation of eqs. (11) and (12), u and v were entered as the mean of 
their depth integral over each of three vertical layers, where the 
depth, Z, of a layer varied across the shelf as a function of the 
bottom depth, H, i.e., Z = H/3. 

2. Light and nutrient regulation 

The spatio-temporal distributions of nitrate (N) and chlorophyll 
(M) were calculated from 
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where u and v are the horizontal velocities, w is the vertical 
velocity derived from the continuity equation, K z is the vertical eddy 
or mixing coefficient discussed below, b is the nitrate/ chlorophyll 
ratio of 0.5, e is the phytoplankton growth rate described in this 
section, w s is the sinking rate of the phytoplankton, and g is the 
grazing loss rate to herbivores. 

In this simulation analysis, we were concerned with the amount of 
phytoplankton export from the mid-Atlantic shelf, which depends on the 
amount of "new" dissolved nitrogen added from the slope and estuarine 
boundaries (EPPLEY and PETERSON, 1980; WALSH, 1981), not the recycled 
production based on ammonium or urea regenerated within the shelf 
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ecosystem. We thus did not consider the growth of phytoplankton in 
eq. (.12) from excretory products (WALSH, 1975), implicit within the 
grazing loss term. 

The input of nitrate from the coastal boundary was assumed to be 
10.0 pg-at NO^ l * from the Hudson River, 1.3 pg-at NO^ £ * from 
Chesapeake Bay, 1.0 pg-at NO^ £ * from Delaware Bay, and 1.0 pg-at NO^ 
£ ^ from Long Island Sound. Most of the model's cases assumed that 
the nitrate boundary condition of slope water and of upstream shelf 
water was 6.0 pg-at NO^ £ * throughout the water column. The 
estuarine boundary conditions of nitrate were later increased tenfold 
for the eutrophication case, while in two simulation runs the upstream 
nitrate content was set instead to 1.0 pg-at NO^ £ 

The specific growth rate (hr *) of the phytoplankton, e, was a 
function of ambient light (I) and nitrate (N) 


e = 


(1 - I/I ) 

e [I/I e s ][ — 5 -jj] 

ms n + N 


(13) 


where e was the maximum growth rate of 0.026 hr , assuming a mean 
m 

temperature of 5°C and a phytoplankton assemblage of 50% netplankton 

(MALONE, 1982). The saturation light intensity, I g , of STEELE's 

-2 -1 

(1962) inhibition expression was taken to be 5 cal cm hr (MALONE, 
1977). The Michaelis half -saturation constant, n, of the nutrient 
limitation term (CAPERON, 1967), at which e = 2, was taken to be 

1 pg-at NO^ £ ^ for coastal phytoplankton species (CARPENTER and 
GUILLARD, 1971). The depth integral of ambient light over each of the 
3 layers was actually a time -dependent function of season, day length, 
diel periodicity, and algal biomass. 
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Such a complex light field was described by 


- -(k + k )z - k Mdz 

l(z,D,t) = 61 sin [ir(t - a)/i(D)]e W s ^ ■'o (14) 

m 

where the seasonal light variation in the model was first described 

(NAGLE, 1978) as a cosine function of incident radiation (I ) at the 

o 

winter solstice, 1^ (120 g cal cm ^ day ^), at the summer solstice, I s 
_2 

(500 g cal cm day ), and the Julian date, D, of the 58-day 
simulation, 

I (D) = I W + 0.5(I S - I W )[1 - cos 2ir(D - 356)/365] (15) 

o o o o 

The changing photoperiod, or day length, was similarly described by 

i(D) = i W + 0.5(i S - i W )[l - cos 2tr(D - 356)/365) (16) 

where i W = 9.1 hr and i S = 15.6 hr (NAGLE, 1978). The 
photosynthetically active radiation (PAR) at depth is about 50%, of 
which 15% might be lost as a result of albedo, such that 6 was taken 
to be 0.425 in this model. 

With specification of the day length, input of total radiation to 
the sea surface, and % transmission of 400-720 nm wavelengths, the 
hourly change of solar irradiance with time, t, can be expressed 
(IKUSHIMA, 1967) as 
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( 17 ) 


I (t) = I sin 3 [n(t - a)/i(D)] 
o m 

where I is the maximum light intensity at local noon (g cal cm 
hr *) for that Julian day, D, and a is the time of sunrise. The solar 
irradiance for the day is 

i(D) 

I (D) = / I sin J [ir(t - a)/i (D) J dt (18) 

o m 

o 


and was used to obtain 1^ in eq. (17), after solution of eqs. (15) and 
(16) for each simulated day of the model. 

Finally, the integrated light field for algal photosynthesis over 
the depth interval, Z, was computed from the self-shading, exponential 
term of eq. (14), where the total extinction coefficient has been 


partitioned into those for water, k w ; detritus, k g ; and phytoplankton. 


kp. The specific attenuation coefficient for plant pigments was 
assumed to be 0.020 m 2 mg" 1 chi a (BANNISTER, 1974; JAMART, WINTER, 


BANSE, ANDERSON and LAM, 1977). The water and detrital contribution 
to water clarity, k w + k s = 0. 13 m 1 , was derived as the average 
residual from the observed extinction coefficient and pigment 
concentration during 15-24 March 1979 (WALSH et al. , 1987a). 

Initial conditions of the spatial chlorophyll field of eqs. (11), 


(12), and (14) were obtained from a CZCS image of the mid-Atlantic 


Bight on 28 February 1979 (SYSTEMS AND APPLIED SCIENCES CORPORATION, 
1984), at about the same spatial resolution as the model's grid 
(Fig. 4). The 1.0 pg 2. 1 surface isopleth of chlorophyll was located 
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near the 20 m isobath, and the 0.5 yg l ^ isopleth at the 60 m 
isobath, with the assumption of uniform pigment distribution within 
the 3 vertical layers of the model. The estuarine boundary conditions 
of algal biomass were respectively 7.5, 5.0, 4.0, and 3.0 yg chi l ^ 
from the Hudson River, Delaware Bay, Long Island Sound, and Chesapeake 
Bay. The chlorophyll contents of inflowing shelf and slope waters 
were instead set to 0.5 and 0.2 yg chi & ^ at the upstream and 
offshore boundaries. 

Since Z was a function of the bottom depth, the upper layer of 
the model was 10 m deep at the 30 m isobath and 30 m deep at the 90 m 
isobath as a result of entering the known bottom topography, H, every 
~3 km in the model. The observed depth of the euphotic zone in the 
mid-Atlantic Bight during March ranges, however, from about 10 m in 
the Hudson River plume (MALONE, HOPKINS, FALKOWSKI and WHITLEDGE, 

1983) to ~30 m at the shelf-break (WALSH et al. , 1987a). Thus a 
surface layer of variable depth (as well as the middle and bottom 
layers) in the model mimics fairly well the spatial changes of 
biological processes across the shelf. 

3. Vertical mixing 

Time series of CZCS images and vertical profiles of chlorophyll 

in March-April 1979 and 1984 (WALSH et al. , 1987a; WALSH et al., 

1987b) implied a downward displacement and/or sinking rate of 20-40 m 

day ^ for biogenic particles. Using a one-dimensional model (NIILER, 

1975) of the wind-induced mixed layer, WROBLEWSKI and RICHMAN (1986) 

2-1 2 

computed a vertical eddy coefficient, K , of 68 ra hr (188 cm 

z 
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sec ^ ) over a 44 m deep mixed layer of weak vertical stratification 
(0.35 o t 50 m *), after 8.5 hr of a 10 m sec ^ wind forcing. The 

2-i 2 -1 

decay time for a return of K to 8 m hr (22 cm sec ) was of the 

z 

same order of 24 hr after cessation of the wind impulse. 

During February-May 1979-82, a wind event £10 m sec ^ occurred 

about every 8 days in the mid-Atlantic Bight and at 27 shelf stations, 

during 16-24 March 1979, the vertical density gradient was a mean of 

0.33 52 m Over a 44 m surface mixed layer in the mid-Atlantic 

2 

Bight, the equivalent vertical displacement rate from a K of 68 i 

z 

hr * would be 37.1 m day ^ during such March wind events. We wished 

to distinguish, however, between downwelling, i.e., w of eqs. (11) and 

(12), vertical mixing, i.e., K , and sinking of phytoplankton, i.e., 

z 

w . 
s 

We thus computed (CSANADY, 1976) smaller values of K from 

z 


-1 

= I£ 

z 200 f 


(19) 


for each of the wind cases (Fig. 5), where x was the mean wind stress 

and p was the density of water. Values of K ranged in this model 

z 

from 17 to 80 cm^ sec i.e., at most about half of WROBLEWSKI and 

RICHMAN'S (1986) estimate, except for one run when all values of K z 

were doubled. Each of the 13 values of K were entered both in 

z 

eqs. (9) and (10) to compute u and v, and in eqs. (11) and (12) to 
compute vertical fluxes of nitrate arid phytoplankton. 
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4. Sinking 

Although laboratory sinking rates of diatoms range from only ~1 
to 10 m day ^ (SMAYDA, 1970), repeated daily field observations of the 
1975 diatom spring bloom at the 80 m isobath in the Baltic Sea 
indicated apparent sinking rates of 30 to 50 m day ^ (BODUNGEN, 
BROCKEL, SMETACEK and ZEITZSCHEL, 1981). Time series of stations 
taken ~2 days apart during the 1978-81 spring blooms in the 
southeastern Bering Sea suggested net iri situ sinking rates at the 
75 m isobath of only ~3-4 m day ^ (WALSH, 1983), but winds were 
favorable for upwelling about 50% of the time. Vertical chlorophyll 
profiles taken every 3 hours at the 70 m isobath off Long Island after 
a 15-20 m sec ^ wind event on 3-6 April 1975 again implied a range in 
sinking velocities of 12-50 m day * (WALSH, 1983)'. Within the present 
model, we employed a range in sinking rates of 0 to 20 m day ^ for w g 
within the water column, at the sea surface, and at the water-sediment 
interface. 


5. Grazing 

Zooplankton grazing rates in the mid-Atlantic Bight during March 
are about 10% of the daily primary production (WALSH, WHITLEDGE, 
BARVENIK, WIRICK, HOWE, ESAIAS and SCOTT, 1978; DAGG and TURNER, 1982; 
SMITH and LANE, 1987), with an increase to ~40% by April (DAGG and 
TURNER, 1982; SMITH and LANE, 1987), and more than 100% by October 
(DAGG and TURNER, 1982). From March to April 1984, the dominant 
herbivore, Calanus f inmarchicus , increased 8-fold in abundance over a 
33-day period (SMITH and LANE, 1987), suggesting an exponential 
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Figure 7. The simulated velocity fields (cm sec" 1 ) of the upper a) 
and lower b) thirds of the water column under a wind 
forcing of 1.07 dynes cm" 2 from 296°T during 5-12 April 



increase of this herbivore population, with a doubling time of about 
10 days. 

In this model, we thus mainly employed an exponential increase of 
the grazing stress, gM, from 28 February to 27 April by using 

g = ln( 1 - G)/[24 - i(D) ] (20) 

where G was 0.03 e ^*^23(59 D) £ or jj > 59 m> while G was a constant 
0.06 for H < 50 m, i.e., there was a spatial gradient of grazing 
pressure. To simulate diurnal migration of the copepod grazers, this 
grazing stress was only imposed on the phytoplankton of the upper and 
middle layers of the model at night. In terms of primary production, 
the exponential grazing stress led to a loss of about 10% of the daily 
growth increment in February and 100% at the end of April, i.e., 
losses to benthos and bacterioplankton were then represented by this 
term as well. In one experiment of this model, no grazing stress was 
assumed, while in another experiment, a linear increase in grazing 
stress from March to April of G = 0.1 + 0.0078(D - 69), for D > 69, 
was instead applied. 

6. Computation 

The circulation sub-model was solved by Gaussian elimination 

(EISENSTAT, SCHULTZ and SHERMAN, 1976), with the steady- state values 

of u, v, w, and K entered in eqs. (11) and (12) for the 13 time 
z 

intervals (Fig. 5). Information about N and M was then obtained, not 
from deduction (CHRISTIE, 1941), but from numerical solution. The 
time -dependent solutions of the biochemical state equations were thus 
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obtained over a staggered grid of the same dimensions of the 
circulation model (Fig. 4), with an upstream or "donor cell" finite 
difference in space (ROACHE, 1976). Integration over time was done 
with a semi- implicit method for the vertical diffusion term and an 
explicit, forward in time, differencing method for the other terms 
(ROACHE, 1976). 


RESULTS 

The simulated flow fields of the circulation model were compared 
(WALSH et al. , 1987a) with AOML current meter observations (MAYER, HAN 
and HANSEN, 1982) for April 1979 at 8 moorings in the New York Bight 
(Fig. 4). During this time period, the barotropic flow was about 90% 
of the transport, agreeing with a prior model's currents within a 
vector error of only ~1 cm sec * and 10° in the New York Bight 
(HOPKINS and DIETERLE, 1987). The extension of this circulation model 
to the mid-Atlantic Bight has been previously discussed (WALSH et al. , 
1987a), and yielded a mean difference between simulated and observed 
currents of 3 cm sec ^ velocity and 35° direction during 5-16 April 
1979. The results of the biochemical calculations will thus be 
stressed in this analysis. 

To provide perspective on the biological response to downwelling 

and upwelling events, however, the flow fields during 1-5 April 

(Fig. 6) and 5-12 April (Fig. 7) are shown for respective wind 

-2 -2 

forcings of 0.63 dynes cm from 063°T and 1.07 dynes cm from 296°T 
(Fig. 5). The surface Ekman layer of the model exhibited onshore 
flows of 10-15 cm sec ^ (Fig. 6a) during the northeast wind forcing of 
1-5 April, for example, while the bottom layer displayed similar 
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offshore flow (Fig. 6b), particularly south of New Jersey. The 

2 -1 

vertical mixing coefficient was computed to be 42.7 cm sec during 
this time period, with maximum downwelling velocities of 23-26 m day * 
on the inner (<30 m depth) and middle (30-60 m depths) parts of the 
shelf. 

With a shift in wind forcing to the northwest during 5-12 April, 

weaker offshore flows of 5-10 cm sec ^ occurred in the surface layer 

of the model (Fig. 7a). Onshore flow of 5-10 cm sec ^ then occurred 

in the lower layer (Fig. 7b), with maximum upwelling velocities of 9 m 

day ^ found on the inner shelf, and as much as 14 m day ^ on the 

2 -1 

middle shelf. A K of 80.9 cm sec , computed during the 5-12 April 

z 

case, at the 44 m isobath would yield an effective mixing velocity of 
-16 m day * to resuspend chlorophyll, if more algal biomass were 
located in the lower layer than in the middle or upper layers of the 
water column. In this situation, a combination of the calculated K 

z 

and w, together with an assumed w g of 20 m day would still allow a 
vertical input of phytoplankton to the euphotic zone at a net rate of 
10 m day ^ . 

Using a constant sinking velocity of 20 m day ^ and the 

time -dependent growth and grazing rates [case (c) of Table 1], we will 

describe the seasonal change in biological response to three such 

downwelling cases on 3 March, 4 April, and 16 April, and three 

upwelling cases on 16 March, 11 April, and 20 April 1979. At the 

beginning of March, for example, the incident radiation in the 

-2 -1 

mid-Atlantic Bight was <250 g cal cm day , which, over a well-mixed 
water column, would yield a mean iji situ light intensity of 
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Figure 10. The CZCS estimate of surface chlorophyll during orbit 2010 
on 18 March 1979. 



respectively 64 and 32 g cal cm day within the surface layers of 

the model at the 30 and 60 m isobaths (Fig. 4). Below an in situ 

-2 -1 

light intensity of 40 g cal cm day , the spring bloom does not 
occur at these latitudes (RILEY, 1967; HITCHCOCK and SMAYDA, 1977). 

Consequently, under a weak northeast wind forcing of 0.21 dynes 
cm ^ from 088°T through 4 March 1979 (Fig. 5), the 0.5 and 1.0 pg 8. ^ 
isopleths of simulated chlorophyll in the surface layer were then 
located along the 20 m isobath (Fig. 8a). The dashed lines of 
Figure 8 and subsequent plots are respectively the 10 m, 20 m, 60 m, 
100 m, and 200 m isobaths. After 4 days of nearshore growth at this 
time in the simulation, the 5.5 pg-at 8. ^ isopleth of nitrate of the 
upper layer was also coincident with the 20 m isobath (Fig. 8b). With 
a sinking velocity of 20 m day and a maximum nearshore ( <30 m 
depth) downwelling velocity of 11.3 m day *,5.0 pg chi 8. * were 
computed in the lower layer near the 20 m isobath (Fig. 8c). The 1 pg 
8. * isopleth of chlorophyll in the lower layer was instead coincident 
with the 60 m isobath, while the 0.5 pg 8. * isopleth was found at the 
100 m isobath, -50-100 km farther seaward than those of the upper 
layer. At a grazing loss of only about 10% of the daily primary 
production, most of this nearshore algal biomass of the model was not 

consumed in early March, but advected downstream and seaward. 

-2 

By the northwest wind forcing of 0.90 dynes cm from 287°T on 
16 March 1979, the 2.5 pg 8, * isopleth of chlorophyll in the lower 
layer had begun to impinge on the 60 m isobath as well (Fig. 9c). 

More than 10 pg chi 8, * was found in the lower layer of the model at 
the 20 m isobath (Fig. 9c), similar to shipboard observations during 
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24 March 1979 (WALSH et al. , 1987a). With this upwelling circulation 

-1 2 -1 
of 10-14 m day and a strong mixing coefficient of 74.2 cm sec , 

the 0.5 yg 1 * isopleth of surface chlorophyll occurred at the 60 m 

isobath (Fig. 9a), similar to the spatial pattern detected on 18 March 

1979 during CZCS orbit 2010 (Fig. 10). After two weeks of growth in 

the model, the nearshore dissolved nitrogen stocks were depleted, with 

the 1 pg-at 2. ^ isopleth of surface nitrate found at the 20 m isobath, 

and the 5.5 pg-at l * isopleth at the 60 m isobath (Fig. 9b). 

Subsequent generations of phytoplankton within shallow waters ( <20 m 

depth) would now grow in the model under nutrient limitation, since 

the ambient nitrate concentrations here were less than the 

half-saturation constant -- recall eq. (13). 

After a month of increasing solar radiation, the iji situ light 

intensity within the upper 30 m at the 90 m isobath was >40 g cal 
-2 

cm in early April, while the nitrate concentration was still 
5.5 pg-at NO^ 2. ^ here (Fig. -lib). The daily grazing stress was now 
about 20% of the outer shelf production, however, and as much as 50% 
of the inner shelf production. Within another downwelling 
circulation on 4 April, i.e., analogous to 3 March, a uniform algal 
biomass of >3 pg chi 2, ^ occurred in surface waters, landward of the 
50 m isobath (Fig. 11a). 

Within the lower layer of the model, which reflects the history 
of the surface primary production, a mid-shelf maximum of >10 pg chi 
2. 1 instead occurred between the 20-60 m isobaths on 4 April 
(Fig. 11c). With less than 0.5 pg-at NO^ 2- * found landward of the 
20 m isobath, except for the estuarine discharge, the dissolved source 
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of particulate nitrogen, sinking out of the surface to nearshore 
bottom waters, had been curtailed by this time in the spring bloom. 

The mid-shelf maximum of the model on 4 April 1979 was also the result 
of prior offshore advection of uneaten phytoplankton from shallow 
waters. A near-bottom maximum of >20 yg chi % ^ was found at 
mid-shelf on the 50-60 m isobaths during a 1-5 April 1984 cruise 
(WALSH etal., 1987b). 

A strong upwelling circulation of 9-14 m day * by 11 April 1979 
(Fig. 7) led to a simulated resuspension of chlorophyll from this 
bottom layer (Fig. 12c) to the top layer (Fig. 12a) at mid-shelf off 
New Jersey to Virginia, and immediately south of Long Island. The 
lower layer of the model was depleted of chlorophyll, while the upper 
layer was enhanced. Thus in situ growth within the euphotic zone 
may not have accounted for all of the surface increment of 
chlorophyll, calculated at mid-shelf between April 4 and 11 
(Fig. 12a). 

Similar mid-shelf maxima of surface algal biomass (Fig. 11a) had 
been detected aboard ship in April 1975 (WALSH et al . , 1978) and 
during orbit 2452 of the CZCS on 19 April 1979 (WALSH et al. , 1987a). 
Although nutrient-rich water was advected shoreward within the aphotic 
zone (Fig. 7) in response to upwelling favorable winds (Fig. 5), the 
nitrate content of the upper layer at mid-shelf was actually reduced 
by ~1 yg-at NO^ £. ^ between April 4 (Fig. lib) and 11 (Fig. 12b). At 
a nitrate/chlorophyll ratio of 0.5, about 40-50% of the primary 
production, i.e., biomass increment + grazing loss, was the result of 
particulate incorporation of dissolved nitrate. 


19 




Another reversal of the winds from the northeast yielded the 

third downwelling case of the model by 16 April 1979, with a K of 

z 

2 -i 

only 31.2 cm sec , but a maximum downward entrainment velocity of 
22-30 m day ^ . High chlorophyll (>10 pg St was again found within 
the bottom layer south of Long Island and at mid-shelf, from New 
Jersey to Delaware (Fig. 13c). In contrast to the second downwelling 
case of 4 April (Fig. 11a), however, a mid-shelf maximum of 2.5 pg chi 
SL * was simulated by 16 April in the surface layer of the model 
(Fig. 13a), with plumes of chlorophyll extending off the Hudson and 
Delaware estuaries. These were also observed by the CZCS during 
orbit 2425 on 17 April 1979 (Fig. 14). Perhaps as a result of an 
apparent lack of surface export of chlorophyll from the Chesapeake 
estuary (Fig. 13a), the chlorophyll content of the model's lower layer 
at mid-shelf off Virginia on 16 April (Fig. 13c) was less than that on 
4 April (Fig. 11c). 

-2 

The third upwelling case, of 0.58 dynes cm on 20 April 1979 

from 341°T, was the weakest northwest wind forcing, with maximum 

positive vertical velocities in the model of only 6-13 m day ^ , and a 
2-1 

K z of 35 cm sec . At the 44 m isobath, such an upwelling rate of 
6 m day ^ and an effective mixing rate of only 7 m day ^ would 
actually result in a net downward algal transfer of 7 m day with 
the assumed sinking rate of 20 m day The simulated surface nitrate 
pattern on 20 April (Fig. 15b) was quite similar to both the previous 
upwelling event on 11 April (Fig. 12b) and the preceding downwelling 
case on 16 April (Fig. 13b), indicating little flux of nitrate to the 
euphotic zone. 
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Figure 14. The CZCS estimate of surface chlorophyll during orbit 2425 
on 17 April 1979. 



A chlorophyll plume of >2.5 pg J l ^ , derived from the Hudson River 
export, was found in the upper layer of the model off New York 
(Fig. 15a), with local maxima of the same amount of algal biomass 
found at mid-shelf, farther to the south. The higher biomass 
simulated previously by the model in the second upwelling case on 
11 April (Fig. 12a) and observed by the CZCS during orbit 2452 on 
19 April 1979 (WALSH et al. , 1987a), was not reproduced by this 
version of the model. Part of the chlorophyll increase detected by 
the CZCS on 19 April may have been the result of surface export of 
macrophytes from coastal marshes, rather than just resuspension of 
near-bottom phytoplankton residues. A lower algal sinking rate of 1 m 
day ^ in the model would have yielded a bloom, of course, in this 
region of the shelf. 

Over the 58-day period of this particular simulation of the 1979 

spring bloom, the phytoplankton export was a mean, over 325 km long 

segments at the shelf -break, of 1.30 g chi m 1 day * between Long 

-2 -1 

Island and New Jersey, 2.25 g chi m day between New Jersey and 

-2 -1 

Maryland, and 3.64 g chi m day between Maryland and Cape Hatteras 
(Fig. 4). More than 90% of this export occurred in the lower layer of 
the model, as indicated by the time history of simulated algal biomass 
at the 80 m isobath, south of Long Island (Fig. 16). This grid point 
of the model was adjacent to 1984 fluorometer moorings at the 80 and 
120 m isobaths (WALSH et al. , 1987b), which exhibited similar records 
of chlorophyll fluctuations at depths of 13 m and 75-81 m (Fig. 17). 
Note that near-bottom chlorophyll concentrations of >2 pg chi t ^ were 
found after Julian day 80 in both the simulated (Fig. 16) and observed 
(Fig. 17) time series. 
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The ability of this case of our simulation model to replicate 

both spatial patterns of chlorophyll detected by the CZCS (Figs. 10, 

14), and temporal patterns detected by moored fluorometers (Fig. 17), 

-2 -1 

suggests that a simulated export of 2.40 g chi m day may be a 

reasonable time-averaged mean for March-April in the mid-Atlantic 

Bight. If we assume that this export occurred at the 120 m isobath, 

with 90% in the lower third of the water column, the depth- averaged 

-2 -1 

mean over the whole water column would then be 0.88 g chi m day 
The possible range and implication of such an estimate is discussed in 
relation to our evaluation of the parameter space of the model. 

DISCUSSION 

Nine separate cases of this simulation model were run (Table 1), 
with, for example, sinking rates of 1 m day * (a), 10 m day * (b), and 
20 m day * (c), at an exponential grazing rate, at the normal mixing 
rate, and at the normal estuarine or upstream boundary conditions. 

The above results were from case (c). In another run, the mixing rate 
was doubled (d), with a sinking rate of 1 m day ^ , and the rest of the 
parameters unchanged. In a third experiment, the zooplankton grazing 
rate was assumed to be zero (e), or to vary linearly with time (f) 
rather than as an exponential function, at an intermediate sinking 
rate of 10 m day with the normal mixing rate and boundary 
conditions. A fourth experiment increased tenfold the estuarine input 
of dissolved nitrogen, at upstream boundary conditions of 6 yg-at NO^ 
l * (g) and 1 yg-at NO^ l * (h), with a sinking rate of 20 m day an 
exponential grazing rate, and the normal mixing rate. Finally, the 
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Table 1. Parameter evaluation of the 1979 spring b loom simulation 
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Figure 16. The simulated chlorophyll content (pg H' 1 ) at noon of each 
Julian day from 1 March to 27 April 1979 within the 
upper middle ( — ), and lower (...) layers of the 

80 m water column at a grid point, adjacent to the 1984 
fluorometer moorings, south of Long Island. 



last experiment considered an upstream boundary condition of only 
1 pg-at NO^ l 1 (i), at the normal estuarine flux of nitrogen, the 
usual mixing rate, a sinking rate of 20 m day ^ , and an exponential 
grazing rate. 

10 2 

With a surface area of the model of 8.52 x 10 m , an upstream 

6 2 . 

shelf boundary area of 8.95 x 10° m , a downstream shelf boundary area 

6 2 7 2 

of 2.10 x 10 m , a shelf -break boundary area of 8.85 x 10 / m , and an 

6 2 

estuarine boundary area of 2.54 x 10° m , the units of Table 1 can be 

g 

converted to total mass fluxes. For example, a mean of 8.9 x 10 g 
chi day ^ was produced over all of the mid-Atlantic shelf during 
March-April 1979 for case (a) with a sinking rate of 1 m day ^ . An 
average of 66% of this synthesized chlorophyll was lost each day to 
grazing, while 21% was lost to export, either downshelf or within 
slope waters. The remaining 13% had not yet been removed from the 
shelf water column by the end of the simulation (Table 1). In 
case (e) of no grazing loss (Table 1), 38% of the March-April daily 
chlorophyll production was instead lost as physical export past the 
shelf and slope boundaries (Fig. 4) over the same time period. 

Using a C/chl ratio of 45/1 (MALONE et al. , 1983), the mean daily 

fixation of carbon, based just on nitrate, would have been 0.47 g C 

_2 -2 -1 
m day in case (a) and 0.31 g C m day in case (c) of the model 

(Table 1). Averaging over upwelling and downwelling wind events, a 

- 2-1 

mean primary production of 0.57 g C m day was actually measured in 
the mid-Atlantic Bight (O'REILLY and BUSCH, 1984) during March-April 
1977-80 at 66 stations, excluding the observations from the mouths of 
the estuaries (WALSH et al. , 1987a). Approximately 50% of the spring 
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CHLOROPHYLL CHLOROPHYLL 



FEBRUARY-APRIL 1984 (JULIAN DAY) 


Figure 17. The observed chlorophyll content (yg f 1 ), after 
application of a 40-hr low-pass filter, from moored 
fluorometers both at a) 13 m and 75 m on the 80 m isobath, 
and at b) 13 m and 81 m on the 120 m isobath, south of 
Long Island from 19 February to 4 April 1984. 


bloom's nitrogen demand is met by nitrate off New York (CONWAY and 

WHITLEDGE, 1979). The total carbon fixation predicted by the model 

-2 -1 -2 
would thus have been 0.94 g C m day in case (a) and 0.62 g C m 

day 1 in case (c), if uptake of recycled nitrogen had been added to 

eq. (12). 

The productivity difference between these two cases of the model 
was the result of increased light limitation, induced by changing the 
sinking rate from 1 m day * to 20 m day ^ (Table 1) i.e., thereby 
decreasing the residence time of an algal cell within the euphotic 
zone. In contrast, a doubling of in case (d) had little impact on 
the production or export of chlorophyll, compared to case (a) in 
Table 1. The vertical structure of chlorophyll in cases (a) or (d) 
was quite different, however, from case (c), either over time at the 
30 m (Fig. 18) and 50 m (Fig. 19) isobaths, south of Delaware, or 
within specific spatial patterns in the downwelling (Fig. 20) and 
upwelling (Fig. 21) circulation modes. 

A sinking rate of 1 m day ^ allowed no vertical structure of 
chlorophyll to develop at the 30 m (Fig. 18a) and 50 m (Fig. 19a) 
isobaths, adjacent to the Philadelphia dump site at 38°37'N, 74°22'W, 
i.e., east of the mouth of Delaware Bay. A uniform 5-6 yg chi l ^ was 
then found throughout the whole mid-shelf region within either the 
lower layer of a downwelling circulation on 4 April 1979 (Fig. 20a)^ or 
in the surface layer during an upwelling event on 11 April 1979 
(Fig. 21a), i.e., Julian days 94-106 of Figure 19a. In cases (a) or 
(d), the seasonal peak of vertically homogeneous chlorophyll on day 95 
at the 50 m isobath (Fig. 19a) followed about two weeks after the 
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60 70 80 90 100 110 


Figure 18. The simulated chlorophyll content (yg i' 1 ) at noon of each 
Julian day from 1 March to 27 April 1979 within the 
upper middle ( — ), and lower (...) layers of the 

30 m water column at a grid point, east of Delaware Bay, 
of a) case (a), b) case (f), and c) case (g). 




maximum algal biomass on day 82 at the 30 m isobath (Fig. 18a), i.e., 

a seaward progression of ~2 km day ^ . 

In cases (a) or (d), a peak of 6 pg chi l * throughout the water 

column was subsequently observed one week later on Julian day 102 at 

the 80 m isobath, south of Long Island, as well. This was in contrast 

to about 4 pg chi l ^ in the lower layer and 1 pg chi l ^ in the upper 

layer of case (c) -- recall Figure 16. Since high values of 6 pg chi 

1 ^ were not found during April from this region, either in the 1984 

moored fluorometer records (Fig. 17), or within a previous 9-yr data 

base of ~5000 stations (Fig. 22), the loss terms of sinking and 

grazing within cases (a) and (d) were probably too small. 

An increased sinking rate of 10 m day ^ and no grazing loss of 

case (e) resulted (Table 1) in the same primary production as 

cases (a) and (d). The average April chlorophyll of the lower layer 

in case (e) was 11.42 pg £ * (Table 2), however, about 3-fold that of 

cases (a) or (d). The average daily algal export to the slope of the 

lower layer of the model past the location of the fluorometer moorings 

-2 -1 

in case (e) also increased to 4.0 g chi m day , compared to an 

observed 1984 mean of 2.7 g chi m ^ day ^ (WALSH et al. , 1987b). The 

computed slope import of case (e) would actually have been 8.0 g chi 
-2 -1 

m day , if recycled nitrogen had been included as a nutrient source 
of this model. 

A linear grazing stress, which was larger than the exponential 
formulation, and the same 10 m day * sinking rate of case (f), instead 
consumed all of the March-April primary production on a mean basis 
(Table 1). The algal biomass above the 30 m (Fig. 18b) and 50 m 
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Figure 19. The simulated chlorophyll content (yg i,' 1 ) at noon of each 
Julian day from 1 March to 27 April 1979 within the 
upper (--), middle ( — ), and lower (...) layers of the 
48 m water column at a grid point, adjacent to the 
Philadelphia dumpsite east of Delaware Bay, of 
a) case (a), b) case (f), and c) case (g). 






(Fig. 19b) isobaths was reduced to zero by Julian day 115 of case (f), 
with only 0. 5-1.0 pig chi l ^ found previously on 4 April (20b) and 
11 April (21b) 1979. The average 1979 March-April export at the 80 m 

isobath of this third grazing experiment in case (f) was only 0.3 g 

-2 -1 

chi m day , an order of magnitude less than the observed flux in 
1984. 

The exponential grazing stress and the 20 m day * sinking rate of 
case (c) thus appeared to be the best combination of loss terms of 
this simulation model to approximate prior satellite (Figs. 10, 14), 
moored fluorometer (Fig. 17), and shipboard (Fig. 21) estimates of 
chlorophyll. We accordingly used these parameters of the model to 
evaluate changes of the boundary fluxes of nutrients within cases (g), 
(h), and (i). The lingering 5 pg-at l ^ concentration of nitrate 
across the shelf, south of Martha's Vineyard, in Figure lib, for 
example, was an artifact of the time invariant upstream boundary 
condition of 6 pg-at NO^ £, ^ in case (c). A reduction of the upstream 
boundary condition to 1 pg-at NO^ l * (Table 1) in cases (h) and (i) 
did not significantly alter the primary production, grazing loss, or 
export, however, compared to case (c), and we discuss instead the 
results of the eutrophication case (g). 

A tenfold increment of the estuarine nutrient loading to the 
mid-Atlantic Bight in case (g) significantly increased the nearshore 
chlorophyll concentrations of the upper (Fig. 21c) and lower 
(Fig. 20c) layers between Cape May and Montauk Point, i.e., the New 
York Bight, as well as off Delaware Bay, and off Chesapeake Bay, 
compared to the previous results of case (c) in Figures 11c and 12a. 


26 





The algal biomass of case (g) did not increase much either at the 30 m 
(Fig. 18c) and 50 m (Fig. 19c) isobaths off Delaware, however, or 
within the whole spatial domain of the model (Table 2). The pr ima ry 
production of the entire mid-Atlantic Bight (Table 1) only increased 
by 5% in case (g). 

Since the export of estuarine nitrogen to the shelf in case (g) 
was still only 10% of that. from the upstream boundary, this result is 
not surprising. It does raise the point, however, that compression of 
4 x 10 7 data points in one single spatio-temporal mean for each entry 
of Table 1 obscures the spatial and temporal complexity of both this 
model and the "real world." Table 1, for example, does not indicate 
for case (c) that either 90% of the simulated slope import occurs in 
the lower layer of the model, or that the computed slope import off 
Virginia is 3-fold that off New York. 

Furthermore, the shelf export has a time dependency associated 
with the estuarine source function. Extension of this calculation for 
case (g) to 20 May 1979, for example, leads to a 25% increase of the 
shelf export from the lower layer, since not all of the 
estuarine -derived particulate matter had reached the shelf-break by 
20 April 1979 in the model. This effect is demonstrated by the first 
and second bar graphs of each month in Figure 21, which are the mean 
chlorophyll concentrations during 1974-82 in the upper 20 m and the 
lower 20-50 m, or 20-75 m, of the New York Bight, while the third and 
fourth bar graphs are other chlorophyll means of an independent data 
set, taken over the larger area of the mid-Atlantic Bight. During 
March, April, and May, the chlorophyll concentrations of the New York 
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Table 2. April 1979 algal biomass (yg chi l *) scenarios 


Experiment 

Upper layer 

Middle layer 

Lower layer 

Water column 

a 

4.01 

3.78 

3.97 

3.92 

b 

2.64 

3.40 

4.87 

3.64 

c 

1.81 

2.90 

5.20 

3.30 

d 

3.95 

3.80 

3.90 

3.88 

e 

6.16 

8.07 

11.42 

8.55 

f 

0.25 

0.33 

0.48 

0.35 

g 

1.97 

3.12 

5.54 

3.55 

h 

1.75 

2.77 

4.85 

3.12 

i 

1.59 

2.54 

4.50 

2.88 





Bight, from coastal to slope waters, are at least twice that of the 
whole mid-Atlantic Bight, reflecting the increased nutrient loading of 
the Hudson River (Fig. 3), i.e., simulated in case (g) as Figures 20c 
and 21c. 

With a decline in nutrient loading from the Hudson River by June 
(Fig. 3) and a seasonal increase of the grazing stress, there are no 
differences in mean algal biomass between the data sets of the New 
York and mid-Atlantic Bights for the months of June to February 
(Fig. 21). If a grazing stress was not imposed during a simulated 
tenfold increase of nutrient export from the Hudson River during 
summer conditions, anoxia, in fact, developed within the New York 
Bight (STODDARD and WALSH, 1987). Temporal and spatial changes of 
biochemical processes within the mid-Atlantic Bight are thus both 
important in determining the fate of primary production within this 
shelf ecosystem. 

The algal biomass left behind in the water column as a net result 
of the birth and death processes of shelf phytoplankton populations 
is, by itself, a poor index of the fate of primary production. 
Averaging over the whole water column, for example, cases (a)-(d) had 
a range in April algal biomass of only 3.30 to 3.92 pg chi il ^ 

(Table 2), but the surface content of case (c) was less than half that 
of case (a). The possible twofold increase of April chlorophyll 
concentrations between 1980 and 1984 (Fig. 1), can thus be simulated 
(Table 2), with a change in the sinking rate of this simple model 
(Table 1). Alternatively, the grazing rate in April 1984 might have 
been less than in April 1980, as indicated by the results of case (e) 
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AVERAGE CHLOROPHYLL (/xg /LITER) 



Figure 22. Monthly distribution of mean chlorophyll (yg i' 1 ) within 
the upper 20 m and the lower 20-50 m, or 20-75 m, of the 
coastal zone, the middle shelf, the outer shelf, and the 
slope in the New York Bight (first and second bar graph of 
each month) and the mid-Atlantic Bight (third and fourth 
graphs) during 1974-1982 from 4673 stations. 













in Table 2. If we had confined the spatial mean of Table 2 to 
estuarine waters ( <20 m) of case (g), moreover, a tenfold increase of 
the nutrient supply would also have led to more than a doubling of the 
algal biomass (Figs. 20c and 21c). 

The three alternative scenarios of light limitation, nutrient 
regulation, and grazing stress are thus each capable of inducing a 
2-fold change in algal biomass within the mid-Atlantic Bight 
(Table 2). In concert, these processes can lead to more than a 
tenfold variation in downstream and offshore shelf export, amounting 
to a range of 8% to 38% of the mean March-April primary production 
(Table 1). We are pleased that such a relatively simple biochemical 
model, coupled to a simple physical model, can both provide insight on 
the interaction of these different processes and reproduce chlorophyll 
fields sensed by satellite, shipboard, and iji situ instruments. 

If recycled nitrogen had been included as a state variable, 

case (c) would have predicted a mean March-April primary production of 
- 2-1 

0.62 g C m day over the whole model domain, with an export of 
-2 -1 

2.60 g chi m day from the lower layer of the model off Long 
Island. Such biological and physical fluxes of organic matter have 
been measured in the mid-Atlantic Bight as well (WALSH et al. , 1987a, 
b). Additional field experiments will provide a robustness to such 
simulation models, which will eventually allow us to predict, rather 
than hindcast, the coastal zone's response to continuing perturbations 
of eutrophication, overfishing, and climate modification. 
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